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DRAFT, February 2, 2008 

Abstract. We have analysed the chemical photospheric composition of 24 Herbig Ae/Be and Vega-type stars in search for the A 
Bootis phenomenon. We present the results of the elemental abundances of the sample stars. Some of the stars were never before 
studied spectroscopically at optical wavelengths. We have determined the projected rotational velocities of our sample stars. 
Furthermore, we discuss stars that depict a (selective) depletion pattern in detail. HD 4881 and HD 139614 seem to display an 
overall deficiency. AB Aur and possibly HD 126367 have subsolar values for the iron abundance, but are almost solar in silicon. 
HD 100546 is the only clear A Bootis star in our sample. 
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1. Introduction 

Herbig Ae/Be (HAEBE) stars are generally thought to be the 
somewhat more massive analogues of the T Tauri stars, which 
are low-mass pre-main-sequence stars. They were first defined 
as a class in 1960 by Herbig, based on three observational cri- 
teria: the stars have spectral type A or earlier, display emission 
lines, and illuminate a bright reflection nebula in their imme- 
diate vicinity. Later on other observational criteria were pro- 
I posed. Finkenzeller & Mundt (1984) mentioned that a shared 
characteristic of Herbig Ae/Be stars is the presence of an IR 
excess, which is an independent, but more easily detected sig- 
nature of circumstellar matter Following The et al. (1994) a 
useful working definition is: Herbig Ae/Be stars have spectral 
type B to F8 and luminosity class III to V, and show emission 
lines. They display a flux excess in the IR, due to the presence of 
cool and/or hot circumstellar matter The observational proper- 
ties of HAEBE stars have been reviewed by e.g. Perez & Grady 
(1997) and Waters & Waelkens (1998). 

That HAEBE stars are young objects has been convincingly 
confirmed by their presence in star-forming regions and also 
by the location in the Hipparcos HR diagram of individual ob- 
jects (e.g. van den Ancker et al. 1998). Some of the objects 
often classified as HAEBE stars may be young main-sequence 
stars rather than pre-main-sequence stars, which is not really 
surprising, since, the more massive the star, the closer the 
circumstellar dissipation timescale approaches the pre-main- 
sequence lifetime. Dusty IR excesses, but without accompa- 
nying gaseous emission lines, are also observed in slightly or 
more evolved main-sequence stars, of which fi Pictoris and 
Vega are the most well known examples. The generic link be- 
tween these debris-disk stars and HAEBE stars is a subject of 



current research to which we aim to contribute with the present 
paper 

Vega is the prototype of the Vega-hke stars and shows mild 
characteristics of the so-called A Bootis stars, which are chem- 
ically peculiar A/F main sequence stars: the elements C, N, 
O and S have approximately solar photospheric abundances, 
while several metals are depleted. Venn & Lambert (1990) 
pointed out that the A Boo abundance pattern is comparable 
to the composition of the gas component of the interstellar 
medium (e.g. Heiter 2002): the depleted elements (Mg, Ca, Ti, 
Fe and Sr) are those with the higher condensation temperatures. 
They suggested that the peculiar abundances could arise from 
gas/dust separation in the circumstellar environment, followed 
by accretion of the metal depleted gas. Reaccretion of circum- 
stellar gas which is 'cleaned' by this mechanism explains the 
most peculiar chemical composition of several post-AGB stars, 
for which the depletion pattern is much more pronounced than 
in the A Bootis stars (e.g. Lambert et al. 1988; Waters et al. 
1992). These peculiar post-AGB stars all turn out to be mem- 
bers of binaries, and the circumstellar dust occurs in a long- 
lived disk which surrounds the binary system (e.g. Van Winckel 
et al. 1995; Waters et al. 1998; Molster et al. 1999). 

Observational evidence (e.g. Mannings & Sargent 1997; 
Testi et al. 2003; Pietu et al. 2003; Fuente et al. 2003; Natta 
et al. 2004) shows that the circumstellar matter in the vicinity 
of HAEBE stars displays a disk-like geometry, the modelling 
of which currently is the subject of intense research (e.g. disks 
powered by accretion (Lin & Papaloizou 1980; Bell & Lin 
1994), irradiated passive disks (Kenyon & Hartmann 1987), ir- 
radiated flared disks (Chiang & Goldreich 1997), disks with 
an inner hole and a puffed up inner rim (Dullemond et al. 
2001; Duflemond 2002; Dullemond & Dominik 2004)). From 
the observation of redshifted resonance lines in the spectra of 
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Table 1. The program stars, listed according to their HD numbers. 
The effective temperature and gravity of the Kurucz model used in 
the present abundance analysis. The last coluirm indicates the cate- 
gory to which the objects belong. References: a-Backman & Paresce 
(1993); b-Bhatt & Manoj (2000); c-Coulson et al. (1998); d-Kalas 
et al. (2002); e-Malfait et al. (1998a); f-Patten & Willson (1991); g- 
Song et al. (2000); h-The et al. (1994). 
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HAEBE stars, evidence for matter accretion has been found. 
These stars then seem natural candidates for displaying the X 
Bootis phenomenon. With the aim of further exploring the pos- 
sible evolutionary connection between the Herbig Ae/Be stars 
and the Vega-type stars, we have carried out a chemical analy- 
sis of a selection of 24 objects, consisting of confirmed HAEBE 
stars and debris-disk stars and other objects which may be of 
interest in this context. 

The plan of the paper is as follows: in the next section we 
describe the sample and the observations; Section 3 is devoted 
to the general results from the chemical analysis and a spe- 
cific discussion of the objects for which depletion may occur. 
In Section 4 we discuss the results. In the Appendix, some notes 
on the individual objects are summarized, as well as a compar- 
ison with Uterature results. 

2. Observations and Analysis 

We have obtained spectra for the 24 objects which are listed 
in Table 1. Fourteen program stars are HAEBE stars, eight 
objects are confirmed Vega- type stars or occur in lists of 
Vega-type candidates, and two stars (HD 18256 and HD 
97633) are rather normal main-sequence A-type stars which 



we measured to check our analysis with respect to litera- 
ture results. A first set of observations, concerning mostly the 
northern-hemisphere objects, was carried out by drs. Hans Van 
Winckel and Gwendolyn Meeus during the nights of 23 and 
25 December 1996, with the William Herschel Telescope and 
the Utrecht Echelle Spectrograph at La Palma Observatory 
(Canary Islands, Spain). The second observational run, con- 
cerning mostly southern- hemisphere objects, was performed 
by dr. Gwendolyn Meeus from January 26 till January 28, 
1999 at the ESQ site La Silla (Chile). These data were col- 
lected using the 1.5m telescope and the FEROS instrument. 
The wavelength range of the spectra lies between 3860 and 
9110A. The spectral resolution of both instruments is compa- 
rable (^ ~ 45, 000). We reduced the data with the MIDAS 
software package. 

Despite the broad wavelength coverage, for several pro- 
gram stars only a few lines turned out to be appropriate for 
detailed abundance analyses. Such a circumstance is not un- 
usual for A/B-type stars, where the condition that spectral lines 
with equivalent widths higher than 150 mA should be avoided 
in abundance analyses and the fairly high projected rotational 
velocities severely limit the amount of suitable lines. For a few 
slowly rotating stars (e.g. HD 18256, HD 20010, HD 104327 
and HD 139614) we were able to use hundreds of lines in the 
analysis; the average value for the amount of lines is about 100. 
It has to be noted that the vast majority of usable lines are iron 
and titanium lines; our search for metal depletion primarily re- 
lies on these two elements, also on silicon and calcium for most 
objects, but on magnesium only in a minority of cases. On the 
other hand, for 6 program stars, no carbon, nitrogen and oxygen 
(CNO) abundances could be determined; the relative scarcity 
of CNO lines and the non-LTE (local thermodynamic equilib- 
rium) sensitivity of especially nitrogen impUes that the uncer- 
tainties on the CNO abundances are sometimes critical for as- 
sessing whether selective depletion occurs; in all cases at most 
a few lines are available for the determination of the sulphur 
abundance. 

As starting values for the stellar parameters (effective tem- 
perature Te and gravity \ogg), we have used previously pub- 
lished results or values derived from photometric data. The 
equivalent widths derived from the spectra were converted 
to abundances using the program MOOG (Sneden 1973) for 
cooler stars (T^e < 8500K) and the Kurucz program W1DTH6 
for hotter stars {T^ > 8500K), both programs supposing LTE. 
Non-LTE effects have not been taken into account. We have 
computed abundances with both programs for the stars with ef- 
fective temperatures in the range from 7500K and 9500K, and 
found good agreement. The stellar parameters were refined - 
when necessary and/or possible- in the usual way: the effective 
temperature was determined by asking that the abundances be 
independent from the excitation potential for ions with many 
observed lines; the gravity was set so that analyses for differ- 
ent ionisation stages for iron led to the same result. The stellar 
parameters as we adopted them in our abundance analyses are 
listed in Table 1 . We were limited in chosing the effective tem- 
perature and gravity parameters by the finite grid of Kurucz 
models at our disposal — the grid has a step of 250K in 
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Table 2. The solar abundances for C, N, O and S are respectively 8.56, 8.05, 8.93 and 7.21 dex (Anders & Grevesse 1989). The given error 
represents the line-to-line scatter. The number of lines that were used in the abundance determination is mentioned in parentheses. 
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and a step of 0.5 dex in log^. In each model, we have used a 
turbulent velocity of Vjur = 3 km s~^ 

3. Results of the chemical analysis 

3.1. Overview 

In this section we present the derived elemental abundances. 
All values are given in reference to the (adopted) solar abun- 
dances given by Anders & Grevesse (1989), [X] - +0.0 in- 
dicating a solar abundance for element X. Since we are in- 
vestigating whether the selective depletion pattern is present 
in our sample stars, we compare the computed abundances 
with the elemental abundances of A Boo and Vega, two pro- 
totype A Bootis stars, and the interstellar medium (Paunzen 
et al. 1999; Venn & Lambert 1990; Qiu et al. 2001; Erspamer 
& North 2002; Jenkins 1989). The stellar parameters for A Boo 
are re=8650K and logg=4.00 (Venn & Lambert 1990) and for 
Vega re=9430K and log ^=3.95 (Qiu et al. 2001). In Table 2 
we have summarized the measured carbon, nitrogen, oxygen 
and sulphur abundances. Table 3 contains the abundances of the 
metals Mg, Si, Ca, Ti, Fe and Sr. In Table 4, we have tabulated 
the abundances we computed for other elements. In parenthe- 
sis, the number of spectral lines that were used in the abun- 
dance determination are indicated^ AH abundances tabulated 
in Tables 2, 3 and 4 and mentioned in the text have a Une-to- 

' The complete linelist and measured equivalent widths can be 
downloaded at 

http://www.ster.kuleuven.ac.be/~brain/LINELIST/linelist.ascii 



Une scatter smaller than 0.20 dex. The typical standard devia- 
tion is A ~ 0.10 dex. The uncertainty on the stellar parameters 
introduces in general another ~0.10 dex. 

Lines -in several cases only one- of additional elements 
could be detected in the spectra of some program stars. The 
results are as follows. The He abundance of HD 250550 
was determined to be [Hei]=-0.03. For HD 20010 we com- 
puted [Mnn]=-i-0.11 dex, which is in good agreement with 
the derived Mm abundance. For HD 104237 we determined 
[A1ii]=-h0.06 dex and [Cui]=-i-0.26 dex (for a solar abundance 
of 4.21 dex). For heavy elements, the following results were 
obtained: the Lan abundances for HD 17206, HD 100453 and 
HD 104237 are h-2.07, h-1.46 and -h0.12 dex, respectively; the 
Ce abundances for HD 17206, HD 18256, HD 33564 and HD 
139614 amount to +1.26, +0.27, -0.15 and +0.64 dex, respec- 
tively; for HD 17206, we found [Ndn]=-0.05 dex, for HD 
33564 [Ndn]=+0.27 dex, and for HD 104237 [Ndn]=+0.59 
dex; the abundance of samarium in HD 104237 is +0.85 dex; 
the europium abundance of HD 17206 is +0.23 dex. The solar 
abundances of lanthanum, cerium, neodymium, samarium and 
europium we used, are respectively 1.22, 1.55, 1.50, 1.00 and 
0.51 dex. 

As a by-product of our spectral analysis, we have computed 
the projected rotational velocity of each star in our sample. The 
values are Usted in Table 5. The typical error is of the order of 
10%. 



Table 3. The abundances of the 'metals'. The adopted solar abundances for Mg, Si, Ca, Ti, Fe and Sr are respectively 7.58, 7.55, 6.36, 4.99, 7.52 and 2.90 dex (Anders & Grevesse 1989). The 
given error represents the line-to-line scatter. The number of lines that were used in the abundance determination is mentioned in parentheses. 
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Table 4. The abundances of other elements. The adopted solar abundances for Na, Al, Sc, V, Cr, Mn, Ni, Zn, Y, Zr and Ba are respectively 6.33, 6.47, 3.10, 4.00, 5.67, 5.39, 6.25, 4.60, 2.24, 2.60 
and 2.13 dex (Anders & Grevesse 1989). The given error represents the line-to-line scatter. The number of lines that were used in the abundance determination is mentioned in parentheses. Two 

numbers refer to the neutral and first ionization state of the element, respectively. 
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HD 31648 


















+0.75 (1) 






HD 33564 


+0.60 (1) 




-0.21±0.10(5) 


-0.61 (1/0) 


-0.21+0.12(5/9) 


-0.29±0.06 (6) 


-0.13±0.11 (11/2) 


-0.56+0.03 (2) 


+0.05+0.11 (3) 


+0.16+0.07(5) 


+0.69 (1) 


HD 36112 






-0.13(1) 




-0.11+0.10(1/6) 


-0.16±0.15 (2) 


-0.18(1/0) 




+0.55+0.11 (2) 






HD 95418 










+0.49+0.06 (0/4) 














HD 97633 


+0.20±0.04 (2) 








+0.17+0.14 (0/7) 














HD 100453 


+0.14+0.14(3) 




-0.21+0.03 (3) 


+0.28 (0/1) 


+0.03+0.10 (3/6) 


-0.04+0.17 (2) 


-0.03+0.15 (7/0) 


-0.52 (1) 




+0.48+0.05 (2) 


-0.02 (1) 


HD 100546 










-0.87+0.05 (0/2) 














HD 104237 


+0.07±0.10(4) 


-0.02+0.12 (4) 


+0.08±0.05 (9) 


-0.01+0.14(1/2) 


+0.17+0.13(16/24) 


-0.05+0.09 (7) 


+0.09+0.14 (34/0) 


-0.37+0.10(2) 


+0.18+0.10(8) 


+0.35+0.11 (6) 


+0.31+0.05 (2) 


HD 135344 






-0.36 (1) 




+0.06±0.06 (2/2) 




+0.08±0.19(2/0) 


-0.45 (1) 








HD 139614 






-0.21±0.08 (5) 


-0.01±0.08(0/2) 


-0.20±0.08 (11/19) 


-0.30±0.01 (2) 


-0.18±0.16(13/0) 


-0.52+0.14(2) 


+0.04+0.18(3) 


+0.47+0.04 (4) 


+0.06+0.01 (2) 


HD 190073 




-0.45 (1) 


-0.05+0.04 (3) 


+0.09+0.08 (3/0) 


+0.01+0.12(2/17) 




+0.06±0.07 (0/2) 




+0.03 (1) 


+0.31 (1) 


+0.11 (1) 


HD 250550 










-0.05 (0/1) 




+0.06 (0/1) 










Vega 
iBoo 
ISM 


+0.12 
-1.30 
-0.50 




-0.77 
-0.70 


-0.26 


-0.48 
-1.00 


-0.30 


+0.02 
-1.20 




-0.59 


-0.93 


-1.22 
-0.90 
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3.2. Stars with evidence for (selective) depletion 
3.2.1. HD4881 

Lagrange-Henri et al. (1990) considered this star, which has 
a Vega-like IR excess, in a list of possible analogues to yS 
Pictoris in a search for redshifted Can and Nai absorption, 
which would point to accretion of circumstellar matter; it is 
not clear from their paper whether they observed the star. 
Miroshnichenko et al. (1999) investigated the optical and IR 
photometry, the Ha emission and measured equivalent widths 
of prominent spectral lines; they found that the star has a far 
IR excess, but rejected it as a Herbig Ae/Be star and classi- 
fied it as a classical Be star. Kalas et al. (2002) discovered a 
Pleiades-like nebulosity around HD 4881, but concluded that 
there was possibly stiU need for a second IR source to explain 
the IRAS fluxes; they suggested that this extra source consists 
of warm grains, either in a circumstellar disk or in a broader 
region where nonequiUbrium small-grain heating occurs. 

We observed HD 4881 in the spectral region from 3860 to 
5020A and from 5220 to 9100A. A few strong spectral fines - 
Mgn A4481, Sin A/14128, 4130, 6347, 6371 and Can /13933- 
in the our spectrum have a peculiar shape (see fig. 1). The pro- 
file seems to be made up of two components — a superposition 
of a broad line and a narrower Gaussian fine. The straightfor- 
ward interpretation is that the narrow line is of circumstellar 
origin, while the broad line is the pressure broadened photo- 
spheric fine. A strange coincidence is then that the broadness 
of the photospheric lines of the other elements in this deficient 
star is close to that of the narrow components of the strong 
fines. The projected rotational velocity we found for HD 4881 
is 25 km s"'. Abt & Morrell (1995) suggested a vsin; of 65 
km s in rather strong disagreement with our value. However, 
Abt & Morrell determined v sin i by fitting a Gaussian curve to 
the /i448 1 Mg n fine, of which the profile is more complex than 
Gaussian. 

Our effective temperature for the Kurucz model of HD 
4881, llOOOK, is consistent with the B - V color index of 
-0.07 (Kalas et al. 2002; Flower 1996) and specfi-al classifi- 
cation B9.5V. Higher effective temperatures lead to a discrep- 
ancy between the Fei and the Fen abundances. We thus dis- 
agree with Miroshnichenko et al. (1999) who proposed a spec- 
tral type B8 for HD 4881 and with the efl'ective temperature of 
12300K derived by Kalas et al. (2002) from fitting the spectral 
energy distribution. Miroshnichenko et al. also determined the 
eff'ective temperature on the basis of the Balmer jump and the 
Ha line, respectively r,=10700±500K and llOOOK; in agree- 
ment with our determination. 

The result of our abundance analysis of HD 4881 is that 
it is deficient in all elements we studied, i.e. this star reveals 
a complete depletion pattern. The iron abundance, measured 
from 4 Fe i fines and 19 Fe n fines, is about -0.70 dex. Similar 
deficiencies are found for titanium and chromium, from 12 and 
9 lines respectively, and for the a element magnesium from 3 
fines, while silicon (5 lines) appears still more deficient. With 
respect to the A Bootis phenomenon, however, we note that 
also for carbon, nitrogen and oxygen subsolar abundances were 
found. 




4470 4480 4490 4500 

wavelength [A] 

Fig. 1. The peculiar shape of the Mg ii line near 4481 A in the spectrum 
of HD 4881. Four Sin lines (4128, 4130, 6347, 6371A) and a Can 
line (3933A) show a similar profile. 



Table 5. The projected rotational velocity. 



Projected rotational velocity 



Star 


V sin i 


Star 


V sin i 


HD 


[kms-'] 


HD 


[km s-'] 


HD4881 


25 


HD 97048 


140 


HD 6028 


95 


HD 97633 


20 


HD 17081 


18 


HD 100453 


39 


HD 17206 


22 


HD 100546 


55 


HD 18256 


17 


HD 102647 


120 


HD 20010 


10 


HD 104237 


10 


HD 28978 


23 


HD 135344 


69 


HD 31293 


100 


HD 139614 


24 


HD 31648 


90 


HD 190073 


9 


HD 33564 


12 


HD 244604 


90 


HD 36112 


60 


HD 250550 


54 


HD 95418 


40 







An alternative explanation for the peculiar line shape ob- 
served in a number of -intrinsically strong- lines may be bina- 
rity. When both components have comparable effective tem- 
peratures, the spectrum of HD 4881 is the superposition of 
two spectra. Under this hypothesis, the line profile displayed 
in Fig. 1 contains a broad feature belonging to the compo- 
nent of the binary with the largest projected rotational veloc- 
ity (vsin J «s 100 km s"'). The narrow absorption feature is 
then due to the second component, which also causes the other 
narrow photospheric absorption lines in the composite spec- 
trum. The binary hypothesis potentially provides a solution for 
the current inconsistencies in determinations of the effective 
temperature and the projected rotational velocity of HD 4881. 
Furthermore, the observed overall chemical depletion pattern 
may be due to veiling effects. Although this hypothesis appears 
very attractive, no independent observational evidence for the 
binary nature of HD 4881 is available in the Uterature. 
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3.2.2. AB Aur 

AB Aur (HD 31293) has been previously mentioned as a A 
Bootis candidate by Renson et al. (1990). Indeed, while we find 
slightly suprasolar values for the nitrogen and oxygen abun- 
dances, iron appears to be strongly depleted. Admittedly, the 
iron abundance is based on only 2 lines, but these lines pro- 
vided reliable abundances for the other sample stars, and even 
taking into account the projected rotational velocity of 100 
km s"\ it seems impossible to assign a solar iron abundance 
to this frequently studied HAEBE star. On the other hand, AB 
Aur appears to be only mildly deficient in sihcon ([Si]=-0.11 
dex). The possible A Bootis nature of this object thus needs 
further investigation. 

3.2.3. HD 100546 

HD 100546 is one of the closest and brightest HAEBE stars. In 
lUE observations, Grady et al. (1995) detected redshifted ab- 
sorption features from iron resonance lines reminiscent of those 
observed for Pictoris and concluded that also HD 100546 un- 
dergoes a bombardment by bodies of cometary size. The ISO 
(Infrared Space Observatory Kessler et al. 1996) spectrum of 
HD 100546 is spectacular, in the sense that it presents promi- 
nent crystalline silicate features with a spectrum closely re- 
sembling that of comet Hale-Bopp (Malfait et al. 1998b). The 
stellar parameters that were used in our Kurucz model for HD 
100546 are taken from Meeus et al. (2001). The projected rota- 
tional velocity of HD 100546 is 55 km s so that most lines in 
the optical spectrum should be suitable for abundance analysis 
purposes. However, it turns out that only a Umited number of 
-intrinsically strong- lines could be measured accurately, and 
that the results consistently point to this star being deficient in 
several important elements. In fact, HD 100546 turns out to 
be the only clear A Bootis star in our sample, being deficient 
about 1 order of magnitude in magnesium (4 lines), silicon (6 
Unes) and iron (5 Unes), but apparently having solar CNO abun- 
dances. Late B stars showing relatively few suitable CNO lines, 
the latter statement is essentially based on our determination of 
the oxygen abundance, which was estimated based on the mea- 
sured equivalent widths of 10 spectral lines -amongst others 
six Unes from the multiplets around A6\56 and /16454 - which 
led to a consistent result with a scatter of 0.10 dex. The nitro- 
gen abundance was computed using only the 7468A Ni line. 
The six Si n lines included the easily recognizable 6347A and 
6371A Unes; the five iron Unes are those at 4233, 4555, 4583, 
5169 and 5275 A, the atomic parameters of which are quite reli- 
able. A similar deficiency was found for chromium, where the 
Unes at /14558 and /14588 were used to compute [Crn]=-0.87 
dex. 



3.2.4. HD 139614 

The star HD 139614 is one of the objects in the sample for 
which the largest amount of Unes can be studied, thanks to a 
small projected rotational velocity, which we found to be 24 
km s~\ in agreement with Dunkin et al. (1997). According to 
Meeus et al. (1998), it is likely that the low v&ini is at least 



A AB Aur 

■ HD 100546 

O A Boo 




5 10 15 20 25 30 

atomic number 

Fig. 2. The elemental abundances of AB Aur and HD 100546 deter- 
mined in the present analysis compared to those of A Boo (Venn & 
Lambert 1990). The depletion pattern of HD 100546 and A Boo are 
clearly similar. This is less obvious for AB Aur. In the lower left cor- 
ner, the typical ±0.20 dex error bar is shown. 



partly due to a near-pole-on incUnation of the star. The stel- 
lar parameters for HD 139614 are again extracted from Meeus 
etal. (2001). 

We found that all elements are slightly deficient ( 0.20 

dex). The abundance pattern of this star is comparable to the 
pattern in HD 4881, but the deficiency is less strong. Dunkin 
et al. (1997) also examined HD 139614 spectroscopically. 
They used a higher effective temperature re=8250K and a 
lower gravity logg=4.2 in their model. They concluded that 
[C]=-0.07 dex, [Mg]=-0.08 dex, [Si]=-0.52 dex, [S]=-0.12 
dex, [Ca]=-0.33 dex, [Ti]=-0.21 dex, [Cr]=-0.11 dex, 
LFeJ=-0.03 dex and [Ni]=-0.07 dex. The abundances of 
carbon, sulphur, calcium, titanium, chromium and nickel are 
in relatively good agreement with our values. The computed 
abundances of magnesium, silicon and iron differ from our 
results by at least 0.15 dex. We note that Dunkin et al. base 
their investigation mostly on strong lines, while we systemat- 
ically discarded all Unes with an equivalent width exceeding 
150 mA. In fact, of the respectively 116 Fei lines and 33 Fen 
lines we used in our analysis, only 1h-4 occur in the Ust of 
lines studied by Dunkin et al. (1997). It is not unUkely that 
our analysis is less afi'ected by the uncertainties induced by 
microturbulence, so that we are confident that the small overall 
underabundance we found for HD 139614 is genuine. 

In Fig. 2 the depletion pattern of AB Aur and HD 100546 is 
compared to the elemental-abundance pattern of A Boo. Fig. 3 
displays the overall-underabundance pattern of HD 448 1 and 
HD 139614. No similarities with the elemental abundances of 
A Boo are noted in these sources. 

4. Discussion 

A few possible mechanisms have been proposed to explain the 
A Bootis abundance pattern. The selective-accretion theory of 
Verm & Lambert (1990) is currently the most widely supported 
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10 20 30 40 

atomic number 

Fig. 3. The elemental abundances of HD 4881 and HD 139614 deter- 
mined in the present analysis. Both patterns display an overall under- 
abundance of -0.65 dex for HD 4881 and -0.22 dex for HD 139614. 
In the lower right comer, the typical +0.20 dex error bar is shown. 

model. In the recent literature, Faraggiana & Bonifacio (1999); 
Gerbaldi et al. (2003) and Faraggiana et al. (2004) mention that 
binarity can play an underestimated role. The composite spec- 
trum might mimic the spectrum of a truly underabundant single 
star. VeiUng due to the companion will make the metallic lines 
in the spectrum appear weaker. It is indeed not unlikely that 
some objects are misclassified as A Bootis stars, since the def- 
inition of this group of objects often used in the literature is 
'A/F stars with weak metallic lines'. We note however that the 
A Boo abundance pattern also requires solar CNO and S abun- 
dances. Veiling due to binarity works as much for the CNO and 
S absorption lines in the spectrum, hence making the measured 
equivalent widths of these lines lower than expected as well. 
The overall abundance pattern however remains roughly unal- 
tered. It is therefore important that metallic as well as CNO and 
S abundances are determined to classify an object as a A Bootis 
star. We stress that the binarity of the selectively depleted post- 
AGB stars cannot account for the observed abundances either. 
This can be deduced from the solar abundance CNO and S, and 
zinc in these objects (e.g. van Winckel 2003). 

Even though young stars with circumstellar (debris) disks 
are a priori fit environments for the selective accretion to take 
place (following Venn & Lambert 1990), the typical abundance 
pattern is not present in the vast majority of these objects. Only 
one star in the present sample, HD 100546, appears to be se- 
lectively depleted. Because of the similarities of the abundance 
pattern of HD 100546 and A Boo, we suggest that accretion of 
metal-depleted gas is a plausible explanation for this source. 
Nevertheless, no other HAEBE or Vega-type star in our sam- 
ple displays this behaviour. A few possible reasons can be pro- 
posed. Selective accretion might be a non-trivial process, which 
only occurs in a minority of the sources. One of the most strin- 
gent assumptions for the phenomenon to happen is the decou- 
pling of the gas and dust particles. When closely interacting, 
the infalling, metal-depleted gas will drag the (metal-rich) dust 
along. Both gas and dust will be accreted in this case, leav- 
ing the photospheric abundances of the central star unaltered. 



A second explanation for the lack of A Boo stars in this sam- 
ple may be the timescales involved. The dust particles are more 
easily prevented from falling in, due to the larger radiation pres- 
sure. Nevertheless, also the dust is expected to eventually spi- 
ral in and be accreted. The depletion pattern is only present 
after the gas has been accreted, and before the dust falls in 
onto the stellar photosphere. The perceptibility of the A Bootis 
phenomenon would then depend on the difference in accretion 
timescale between gas and dust particles. Moreover, Turcotte 
& Charbonneau (1993) have modeled the effect of stellar ro- 
tation and chemical diflftision on the photospheric abundance 
pattern. They find that a constant infall of at least 10"'"* M© yr"' 
is needed to maintain the depletion of the metals when photo- 
spheric diffusion is taken into account. After accretion of gas 
has stopped, the abundance pattern is destroyed on a timescale 
of 10^ yr. Stellar rotation and the consequent meridional circu- 
lation only plays a role for rotational velocities above ~ 125 
km s"'. 

Imaging of HD 100546 shows that the inclination of the 
circumstellar disk around this source is i ~ 50° (Pantin et al. 
2000; Augereau et al. 2001; Grady et al. 2001) where / = 0° 
represents pole-on. Assuming that the rotational vector of the 
star and the disk are located along the same axis, the rotational 
velocity of HD 100546 is v 70 km s~'. This is consistent 
with Turcotte & Charbonneau (1993), since the A Bootis abun- 
dance pattern is observable. AB Aur on the other hand is a 
rapidly rotating star. Using the inclination of / < 30° (Eisner 
et al. 2003; Blake & Boogert 2004) of the circumstellar disk 
for the stellar rotation, the rotational velocity for this object is 
V > 200 km s In this star, the photospheric nitrogen, oxygen 
and silicon are not underabundant, indicating a normal abun- 
dance pattern. Strangely enough, iron is clearly depleted, and 
this signature has not been wiped out. Altogether, fast rotation 
alone cannot account for the lack of A Boo stars in this sam- 
ple; assuming an average inclination angle of / = 45°, only 5 
sample sources (22%) are expected to have rotational velocities 
above 125 km s"'. 

HAEBE stars have ages of a few Myr and mass accretion is 
ongoing in a large fraction of these sources. When the dust par- 
ticles in the circumstellar disk coagulate and grow, the coupling 
between gas and dust becomes smaller. Therefore, one intu- 
itively expects that the A Bootis abundance pattern, which only 
occurs under stringent conditions, appears in the latest phase 
of the pre-main-sequence evolution of HAEBE stars (or later). 
During earlier stages, the mass infall might be too vicious, 
dragging along gas and dust. A necessary, but not satisfactory, 
condition for the depletion pattern to occur seems to be a quiet 
circumstellar environment, in which the dust has grown to ap- 
preciable size or disappeared, and a constant low-mass flow of 
depleted gas accretes onto the central star. 

Vega-type stars are believed to be the successors of HAEBE 
stars. The circumstellar matter around these stars consists of the 
left-overs of the pre-main-sequence dust disk (cold particles at 
large distances to the central star and possibly a planetary sys- 
tem). None of our sample stars belonging to this group displays 
selective depletion of the metals. Gas accretion in most of these 
sources may have become too weak to produce a clear under- 
abundance pattern. 
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HD 100546 is probably at the end of its pre-main-sequence 
lifetime. The observed A Bootis-like abundance pattern adds 
another pecularity to the list of curiosa (e.g., high degree of 
silicate crystallinity (Malfait et al. 1998b), spiral structure in 
the circumstellar disk (Grady et al. 2001)) known about this 
source. 

Appendix: notes on individual objects 

HD 6028. Stellar parameters come from Coulson et al. (1998). 
HD 17081. Adehnan et al. (1991) determined the iron 

abundance [Fe]=-0.07 dex. Takeda et al. (1999) determined 
[OJ=-0.14 dex, [Fe]=-0.02 dex and v sin 19. 2 km s ^ Our 
iron abundance is ~0.20 dex lower. The oxygen abundance and 
the projected rotational velocity are in good agreement with 
our values. Smith & Dworetsky (1993) suggested [Cr]=+0.33 
dex, [Fe]=-0.02 dex and [Ni]=-0.05 dex. Only the nickel 
abundance is in agreement with our determinations. 
HD 17206. Marsakov & Shevelev (1995) derived an iron 
abundance of -0.07 dex. This is in excellent agreement with 
our results. 

HD 18256. Eggen (1998) computed [Fe]=+0.13 dex which is 
fairly close to our determination. 

HD 20010. Santos et al. (2001) and Marsakov & Shevelev 
(1995) derived respectively [Fe]=-0.20 dex and [Fe]=-0.23 
dex . This strokes with our results. 

HD 28978. The computed abundances of Lemke (1989, 
1990) [Ci]=-0.05 dex, [Sin]=-0.06 dex, [Cai]=-0.07 dex, 
[Tin]=+0.15 dex, [Fei]=+0.07 dex, [Fen]=+0.19 dex and 
[Ban]=+0.61 dex are in fairly good agreement with our 
results. The result for the heavy element strontium differs 
from our determinations: [Srn]Lemke=+0-55 dex versus our 
[Srn]=+0.31 dex. 

HD 31648. Stellar parameters are taken from van den Ancker 
et al. (1998). Our measured projected velocity agrees reason- 
ably with the value v sin ;=102 km s^'of Mora et al. (2001). 
HD 33564. DeUyannis et al. (1998) suggested [Fe]=-0.05 
dex. Marsakov & Shevelev (1995) derived [Fe]=-0.11 dex. 
We find a stronger underabundance of iron, [Fei]=-0.35 dex 
and [FenJ=-0.41 dex. 

HD 36112. Beskrovnaya et al. (1999) determined [Si]ai+0.40 
dex, [Ca]=+0.06 dex, [Fe]=-0.02 dex and vsin/=60 km s"'. 
The calcium and iron abundances and the projected rotational 
velocity are consistent with our findings. The silicon abun- 
dance of Beskrovnaya et al. is ~0.40 dex higher than our value. 
HD 95418. Hill (1995) derived [Fe]=+0.14 dex which is 
comparable to our results. 

HD 97048. Stellar parameters come from van den Ancker 
et al. (1998). Because of the large projected rotational velocity 
(v sin i=140 km s~') the equivalent widths of the spectral lines 
are difficult to determine. We based our oxygen abundance 
computation on two O i multiplets at 6156A and 6454A. 
HD 97633. Lemke (1989, 1990) determined [Sin]=+0.15 
dex, [Cai]=-H0.17 dex, [Tin]=+0.20 dex, [Fei]=-H0.24 dex 
and [Fen] =+0.30 dex. The silicon abundance is in good 
agreement with our value. The other elemental abundances 
are systematically ~ 0.20 dex higher than our results. Our 
iron abundance ([Fe]=-i-0.02 dex) confirms the result of 



Adelman (1988). Takeda et al. (1999) computed [O]=-0.27 
dex, [Fe]=-i-0.19 dex, [Si]=-0.43 dex and v sin j=21.6 km s'K 
There is a discrepancy with our silicon and iron abundances. 
Smith & Dworetsky (1993) determined [Cr]=+0.43 dex and 
[Fe]=+0.03 dex. We find a chromium abundance of +0.17 dex. 
HD 100453. Stellar parameters were taken from Meeus et al. 
(2001). We altered the log^ value from 4.5 to 4.0. About 180 
spectral lines were used in the abundance determinations. 
HD 102647. Stellar parameters were extracted from Habing 
et al. (2000). 

HD 104237. We took the stellar parameters of Meeus et al. 
(2001). Close to 440 lines were included in the abundance 
computations. 

HD 135344. We adopted the gravity parameter of Meeus et al. 
(2001) from log ^=4.5 down to 4.0. We kept from that 
article. The projected rotational velocity is 80 km s"^ We 
specifically looked for unblended spectral lines by comparing 
the spectrum of HD 135344 with the spectra of HD 18256 
and HD 20010, two stars with a similar efi'ective temperature, 
but much lower v sin i. By doing this we diminish the risk of 
misidentifying lines and measuring wrong equivalent widths. 
HD 190073. Based on the A3111e spectral classification by 
Gray & Corbally (1998) and the A2IVe classification by Mora 
et al. (2001) we chose an efifective temperature T^. = 9250 K 
and a gravity logg = 3.5 for the Kurucz model. Other articles 
suggested a much higher efl'ective temperature (e.g. 11120K 
and a B91V classification, Cidale et al. (2001) or AOIV by 
The et al. (1994)). Nevertheless, the lower temperature was 
consistent with the spectrum of HD 190073 (see fig. 4). Some 
of the spectral lines of Ca i, Ca ii, Sc ii, Ti ii, Crn, Fe i and Fe n 
depict emission features. The emission fills the core of the 
photospheric absorption line. Hence, the equivalent width of 
the line cannot be measured accurately. The spectral lines that 
display emission are therefore not included in the abundance 
analysis ofHD 190073. 

HD 244604. The stellar parameters of HD 244604 are ex- 
tracted from Miroshnichenko et al. (1999). 
HD 245185. Not one single absorption line was detected. 
Our signal-to-noise level in this spectrum was about S/N~30, 
which makes it by far the worst spectrum in our sample. The 
projected rotational velocity could not be derived. Stellar 
parameters for HD 245185 come from Testi et al. (1998). 
HD 250550. Stellar parameters were taken from Testi et al. 
(1998). 
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Fig. 4. The spectra of HD 28978 (top), an A2 star with r,=9250K, and 
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